clear;
clc;
x1=[0:pi/30:pi]';
y=[sin(x1) cos(x1) sin(x1/2)];
S=trapz(x1,y)